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A recently proposed master equation in the Lindblad form is studied with respect to covariance 
properties and existence of a stationary solution. The master equation describes the interaction 
of a test particle with a quantum fluid, the so-called Rayleigh gas, and is characterized by the 
appearance of a two-point correlation function known as dynamic structure factor, which reflects 
symmetry and statistical mechanics properties of the fluid. In the case of a free gas all relevant 
physical parameters, such as fugacity, ratio between the masses, momentum transfer and energy 
transfer are put into evidence, giving an exact expansion of the dynamic structure factor. The limit 
in which these quantities are small is then considered. In particular in the Brownian limit a Fokker- 
Planck equation is obtained in which the corrections due to quantum statistics can be explicitly 
(— **; ■ evaluated and are given in terms of the Bose function go(z) and the Fermi function fo{z s 
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I. INTRODUCTION 



The study of dissipative systems in a quantum mechanical framework is a subject of major interest for many physical 
communities especially in connection with applications. Recently however the subject has gained new interest also 
for physicists concerned with foundations of quantum mechanics, due to the relevance of the notion of decoherence 
f~^ ' as a gateway between classical and quantum world pF. This interest is strongly supported and partially motivated 
by a spectacular improvement in many experimental techniques useful for handling with great precision single- or 
few- particle systems, checking for coherence properties in their dynamical evolution Bj. In this connection models 
for quantum dissipation determined by the symmetry properties of the microphysical interaction and by symmetry 
r"| ■ and statistical mechanics properties of the environment could be of interest for a large class of phenomena. In 
<*S || the Markovian limit quantum-dynamical semigroups [|3| seem the most natural quantum mechanical framework for 
, 'i . the description of dissipative systems B and a lot of work has been done in this direction, both at rigorous and 
formal level, especially with reference to the structural result of Lindblad, which fixes the form of the generator 
of a completely positive quantum-dynamical semigroup in the case of norm continuity B. In this paper we will 
study in detail some structural properties of a recently proposed Markovian master equation for the description of 
the dynamics of a test particle in a fluid P-fJ, the so-called Rayleigh gas B. This simple but realistic model is 
of particular interest in statistical mechanics, being a paradigmatic example which opens the way to the study of 
interacting truly many-body systems. In the quantum regime a feature of independent interest is the relevance of 
quantum statistics of particles making up the fluid, especially in connection with the recent experimental realization of 
almost non-interacting degenerate gas samples of both Bose and Fermi particles (Hi . The considered master equation 
was derived by assuming a translationally invariant interaction between the test particle and a homogeneous fluid, 
and has the general structure of generator of translation-covariant dynamical semigroups considered in |ll| l . In this 
way a direct physical interpretation arises, at least in a particular case, of the general structure given in ]ll[ |. Further 
symmetry or equilibrium properties, which are of fundamental relevance in order to determine the realm of validity of 
a given master equation [pi, are shown to be a direct consequence of particular physical features of the environment, 
embodied in a specific two-point correlation function, the so-called dynamic structure factor |12| . 

The paper is organized as follows: in Section O the master equation is introduced and the property of covariance with 
respect to translations and rotations is considered with reference to the corresponding symmetries of the environment 
embodied in its dynamic structure factor, moreover the existence of a stationary solution is proved, provided the 



environment is in a /3-KMS state |13|; in Section III the case of a free quantum gas is considered, the dynamic structure 



factor is explicitly calculated, its exact expansion with respect to the relevant physical parameters is obtained and in 



the limit of small fugacity z the expression for Maxwell-Boltzmann particles is recovered; in Section IV the Brownian 
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limit in which the test particle is much heavier than the particles making up the gas is dealt with, together with the 
limit of small momentum transfer, leading from the master equation to a Fokker-Planck equation, strongly dependent 
on the statistics of the gas; in Section |y| the obtained results are briefly summarized and discussed. 

II. GENERAL FEATURES OF THE MASTER EQUATION 

Let us recall the general expression of the master equation proposed in MM for the description of the motion of 
a test particle in a homogeneous fluid supposed to be at equilibrium, whose properties we are going to study. The 
obtained result is based on a scattering theory derivation, assuming a translationally invariant interaction in terms 
of two-particle collisions, and is expected to be valid on a time scale much longer than the relaxation time of the 
environment. In the Schrodinger picture the master equation is given by 

§ = M[Q] = -jr[H ,&]+m a) 

where Ho = p /2M is the Hamiltonian of the free particle, M being the mass of the test particle, while the dissipative 
part is given by the following mapping with a Lindblad structure 



JR 3 Z 



(2) 



where the integral is over the parameter space of the translation group in momentum space in three dimensions, the 
parameter q being the momentum transferred in a collision, {A, B} denotes the Jordan product A o B = AB + BA, 
and the operator p is the generator of translations. The unitary operators U(q) = e^ qx , q £ R 3 , are the generators 
of translations in momentum space or boosts. The positive measure dfi(q) is given by 

du(q) = ^(2irh) 3 n\i(q)\ 2 d 3 q, (3) 

thus being invariant under both rotations and translations. In fact n is the particle density in the macroscopic system, 
and the function t(q), given by 

J R 3 {2Trhy 

where q — \q\ and x — \x\, is the Fourier transform with respect to the transferred momentum q of the T matrix 
describing the translationally and rotationally invariant interaction between the test particle and the particles of the 
fluid, which is supposed to be energy independent. The function S(q,p), which appears operator- valued in (|2|), is a 
positive two-point correlation function known as dynamic structure factor pjfl , given by the Fourier transform with 
respect to energy transfer E{q,p) and momentum transfer q of the time dependent pair correlation function of the 
fluid 

S(q,p) = -!-r f dt [ d^xe^^i^-i-^Gix^), (4) 

where G(x, t) is the time dependent pair correlation function 

G(x,t) = ±=[ d 3 y(N(y)N(y + x 7 t)), (5) 

N(y) being the operator density of particles in the fluid and (. . .) the ensemble average with respect to the state of 
the macroscopic system. The expression E(q,p) gives the energy transfer in a collision where the test particle changes 
its momentum from p to p + q, so that 

E(q,p) = E p+q -E p = ^ + ^. (6) 

Note that in (0) we have used as variables momentum and energy transferred to the test particle. In the sequel we 
will use both the equivalent notations S(q,p) and S{q,E{q 1 p)) = S(q,E), according to convenience. The dynamic 



structure factor is a very important physical quantity, giving the spectrum of spontaneous fluctuations of the system, 
and it is of direct experimental access: in fact, as first obtained by van Hove in a fundamental work |14[] , it is directly 
related to the energy dependent differential cross-section per target particle describing scattering of a microscopic 
probe off a macroscopic sample through the formula 

d 2 a . t . 6 / M V p' ... 



referring to scattering of the microscopic probe from p top' = p + q. The dynamic structure factor is given in (Q) 
as a function of momentum and energy transfer, which are the measured quantities in scattering experiments, the 
energy E being related to q and p through (g). The appearance of the dynamic structure factor in (0) gives the 
physical reason for its being positive definite for every system. The main point in (|l|) and (g) is the determination 
of the specific expressions for the measure dfJ,(q) given by (|3[), and of the operator valued function S(q, •), given by 
(||), which can only be obtained on the basis of a microphysical derivation of the equation, relying on some physical 
model. A general structure encompassing (0) has already been considered by Holevo in a purely mathematical context, 
studying the general structure of generators of translation-covariant dynamical semigroups |0] . In particular Holevo 
has proven that when the generator is bounded it must have a structure of the form (uy with C given by (0) provided 
all operators appearing in ([!]) and (0) are bounded [[JT[ (see also |15],[f6[ , where further restrictions to the structure of 
(0) appear, to be discussed later on), and allowing, instead of the structure 



appearing in (|2|), where S(q,p) is self-adjoint and positive, the more general structure 

U(q)L(q,p)-L\q,p)W(q). 

If the generator is unbounded also diffusion terms of the form considered in ( p7|) may appear, and the operators 
appearing in (||) and (||) may be unbounded (see |ll[ for further details). In the general model considered here the 
Hamiltonian is given by the unbounded operator Hq = p /2M, while the remaining part of the generator is determined 
by the explicit expression of the physical quantities t(q) and S(q,p), depending on the specific model for the fluid. 

We now consider the behavior of O) with respect to symmetry transformations. Let us consider a locally-compact 
group G and a unitary representation U(g), with g £ G, on the Hilbert space of the system. Following [pTj we say 
that a mapping M. in the Schrodinger picture is G-covariant if it commutes with the mapping U g [-] — U(g) ■ U'(g) for 
all g e G 

M[U g []}=U g [M[}}. (8) 

Let us first consider the case of spatial translations. Then the unitary operators are given by U(a) = e _ n ap with 
a E R 3 and exploiting 



[U(a),H ]=0, [U(a),S(q,p)}=0, 



together with the Weyl CCR 



U(q)U(a) = U(a)U(q)e* 



aq 



one immediately has that the mapping M. given by (0) is translation-covariant. This property goes back to homo- 
geneity of the fluid and translational invariance of the interaction, as can be seen in the derivation of the master 
equation [Mpl. We then consider invariance under rotations, so that the relevant set of unitary operators takes the 
form U(R), with R £ SO(3). In this case exploiting [U(R),Hq] — and the relations 

U f {R)U(q)U(R) = U(R- l q), U f {R)pU{R) = Rp, 
one has that M. is rotation-covariant provided the dynamic structure factor satisfies for R £ SO(3) 

S(Rq,Rp) = S(q,p). (9) 

In fact if (0) holds one has 



M[U R [-}]= / dn(q)U(R) 

JR 3 



Mq) u(R) 



U(R- 1 q)VS(q,Rp)- i/sfcm&iR-iq) - ~ {S{q, Rp), ■} 



uHr) 



/r 3 
■U R [M[]]. 



U(q)^S(qJ>)-VS(q~p)U\q) - - {S(q,p), ■} 



uHR) 



On the other hand (M) is directly linked to rotational invariance of the surrounding environment, as one can see 
observing that E(q,p) as given by rtq) satisfies 



E(Rq,Rp)=E(q,p), 



and consider the identity 



S(Rq,Rp) = — f dt f (fxe^^i^P^-^^Gix.t) 

= -Ki dtf d 3 xe^ [E{q ' p)t - qx] G{Rx,t) 
^h J R J R3 

so that (0) holds if and only if the pair correlation function is invariant under rotations, 

G(Rx,t) = G(x,t). 

In order to proceed further and consider the existence of stationary solutions we make the natural assumption that 
the state of the macroscopic system, with respect to which the expectation value in (0) is calculated, is a /3-KMS 
state, so that the relation 



(A(w)B) = (BA(w + ih(3)} 



(10) 



holds. This in turn implies that the dynamic structure factor S(q,p) satisfies an identity known as detailed balance 
condition [[121 



S(q,p)=e-K™ + ^S(-q,p + q), 
usually expressed in terms of the dependence on the transferred energy 

S(q,E(q,p)) = e-^^S(-q,-E(q,p)), 



(11) 



(12) 



the sign of the exponential being determined by the fact that we are considering momentum and energy transferred 
to the particle. It will prove useful for further considerations to introduce a symmetrized version of the dynamic 
structure factor, given by 



S(q,E) = e* E S(q,E) 



or equivalently 



1 2Af ^ M J 



S(q,p), 



S(q,p) = e' 
satisfying instead of (|l2j) the more symmetric 

S(q,E)=S(-q,-E), 
so that (O) becomes 

S(q,p) = S(-q,p + q). 



(13) 



(14) 



(15) 



We now look for a stationary solution of (II), given that the environment is in a /3-KMS state, so that due to (fit]) 
the dynamic structure factor satisfies the detailed balance condition as in (O). According to translation-covariance 



of the generators we look for a solution invariant under translation, of the form p(p). Since [Hq, p(p)] = 0, p(p) will 
be a stationary solution of (hi) provided 



C[P(P)} 



R ! 



dp(q) 



R ! 



dp{q) 



U{q)S{q,p)p{p)U\q) - S(q,p)p(p) 
S(q, P - q)p(p -q)~ S(q, p)p{p) 



(16) 



= 0. 



Introducing the function 



the requirement (16) becomes 



A(q,p) = S(q, p - q)p(p - q) - S(q, p)p(p) 



f dp(q)A(q,p) = Q. 

iR 3 



(17) 



(18) 



A sufficient condition for (18) to be valid is that A(q,p) be an odd function in q, and we shall see that this is exactly 
the case if p(p) has the canonical structure po(P) = e _|5m , with M the mass of the test particle and (3 the inverse 

- 2 

temperature of the macroscopic system, as it is to be expected on physical grounds. Let po(p) = e~^™" ( then 



Po{p- q) = Po(p)e 0{m 



so that 



A(q,p) = Po(p) 



S(q,P~q)e ^% M ^) - S(q,p) 



and exploiting flll| ) 



A{q,p) = Po{p) [S(-q,p) - S{q,p)} , 

which is manifestly odd in q. 

As mentioned above the structure given in (Q) is a particular realization of the general expression considered 
in Jp^l, however it does not meet the more stringent requirements exhibited by the dissipative mapping considered 
in |16| . According to mU these requirements are unnecessary if one is looking for the most general translation-covariant 
generator and we shall show that in the present framework they would lead to unphysical results. In fact the structure 
proposed in Jl6| for the dissipative part would take in the Schrodinger picture the form 



C\ 



with the further requirement 



R :1 



dp(q) 



P{q,p) ■ pHqiP) - 2 {P\q,p)P{q,p), ■} 



p f (q,p) = P(-q,P), 



which in our case, since 

P(q,p) = U(q)yfS(q~^p~) 
is not satisfied, being equivalent to the requirement 



(19) 



(20) 



(21) 



U{q)VS{^P~) = U(q)\/S(-q,P + q), 

which does not hold due to the presence of the factor e~ /3 ^2M+^r) in (O). If instead of the dynamic structure factor 
S(q,p) one would consider the symmetrized dynamic structure factor S(q,p) given by (114) , so that P(q,p) in ( |2l| ) 
would be replaced by 



P(q,p) = U(g)sJS(q,p), 

then according to ( Jig ) the relation (EQ) would hold and the dissipative mapping L would conform to the structure 

proposed in jT6|. In this case, however, one would have the unphysical result that po(P) = 
stationary solution, because 

A(q,p) = S(q,p - q)p (p - q) - S(q,p)p (p) 



w is no more a 



is no more odd in q. Other unphysical features linked to the further restriction fl20|) will be considered in Section IV. 
Note that these features can only be discovered with reference to a specific structure of translation-covariant generator 
determined by starting from some microphysical model. In fact the restriction (|20|), though unnecessary from a 
mathematical standpoint, could have proven as well interesting from a physical point of view, thus suggesting the 
result of jlq] as a useful starting point for phenomenological approaches. This does not seem to be the case, in fact 
the substitution S(q,p) — > S(q,p), natural in order to comply with pi, leads to unphysical results. 



III. EXACT EXPRESSION FOR A FREE QUANTUM GAS 

It is of course of interest to analyze the mapping M. given in (fil) for a model in which the dynamic structure factor 
of the fluid can be explicitly calculated: this is the case for an ideal quantum gas considered at finite temperature 
T = l//3k, where k is Boltzmann's constant, and obeying either Bose or Fermi statistics. Apart from simplicity, the 
case of a free gas can be of interest also in view of the recent experimental realization of dilute quantum samples 
of Bose or Fermi particles in the degenerate regime jL0| . The dynamic structure factor for an ideal gas takes the 
form 111 



S B /v(q,p) 



1 



n J R 3 (2nh) 



d3k I \ M-U/ \ uAP + ?) 2 

3 (nfe) B/ p(l ± (nfc_g)B/FJd 



(fe 



2M 



2m 



P 

2M 



2m 



(22) 



where the indexes B or F and signs + or — refer to Bose or Fermi statistics respectively, M is the mass of the test 
particle, m is the mass of the particles making up the gas, n is the density of particles in the gas and (rifc) B/F is 



[Tlk/B/F 



z -l e f3e k zp 1' 



£fc 



Jl. 

2m 



where z is the fugacity of the gas, related to the chemical potential p by z = e^. For a Bose gas at finite temperature 
< z < 1, while f or a Fermi gas z > 0. The integral in (|22| ) can be explicitly calculated both for bosons and fermions, 
giving the result (A6) obtained in Appendix |AJ: 



SWf(<Z,p) =T 



1 2nm 2 



(2irh) 3 nl3q \-efi E ^,p) 



log 



1 =F z exp 


(3 (2m£(q,p)+g 2 ) 2 l 


8m q 2 


1 =F z exp 


(3 (2mE(q,p)-g 2 ) 2 " 


8m q 2 



(23) 



with q = \q\. In the same way one can consider the case of a free gas of particles satisfying Maxwell-Boltzmann 
statistics, thus having 



SMs(q,p) = 



1 



d\ 



n J R 3 (2nh) 3 



(nfe) M B<5 



(p + q) 

2M 



(k 



P 



2m 



k '.\ 



2M 2m ' 



(24) 



with 



so that the integral in ( p4|) can also be explicitly calculated giving the expression ( |Aq) 



S MB (q,p) = 



1 2nm 2 
(2ttT%Y n/3q 



■zexp 



/? (2mE(q,p) + q 2 ) 2 
8m q 2 



(25) 



A convenient way to write (E3J) and (|25|) for later expansions is in terms of the function 



a(q,p) = — [q 2 + 2aME(q,p)] 



(26) 



where the ratio a = m/M between the masses of the particles of the gas and of the test particle has been put into 
evidence, thus obtaining respectively 



S B/F (q,p) = T 



1 2nm 2 



1 



(2nh) s n(3q X _ exp JL (2a(q,p)q - q 2 ) 



T lo S 



1 =p z exp 



£° 2 (*,p) 



1 =p z exp 



-£W<z,p)-?) 2 



and 



S MB (q,p) 



1 2nm 2 
(2ttT%Y n(3q 



■zexp 



P 2 



2m 



(r{q,p) 



(27) 



(28) 



We have thus put into evidence all the physical parameters which are of interest in specifying the physical model 
under consideration and its range of validity: q, E, a and z. In this perspective the expression for a gas of Maxwell- 
Boltzmann particles can also be obtained as expected from the dynamic structure factor for a Bose or Fermi gas in 
the limit of small fugacity z. In fact starting from ( |27| ) and expanding the logarithm up to first order in z one has 

r. i „ 1 2-kvh 2 z 

S B/F {q,p,z<^ 1) = 



(2irh) 3 n(3q 1 _ exp 


'^(Mq,p)<i-<i 2 ) 








x < exp 


'-|V(«,p) 


— exp 


^ II A 


-if 


1 2nm 2 


':t^"\ 


= S MB (q,p). 




{2irhf nf3q " °^ P 





Both (|23[) and (|2J])or equivalently ( |27| ) and (f28|) are invariant under rotation, as one can see from the fact that they 
depend on q and p only through E(q,p) and the modulus q of q, so that 

S B/F (Rq, Rp) = S B/F (q,p), S MB (Rq, Rp) = 5 MB (q, p) , 

thus leading to a rotation-covariant mapping M. when substituted in (pi). In order to grant the existence of the 
stationary solution po(p) we have to check that the obtained expressions satisfy the principle of detailed balance. 
Starting from (|23|) we have, setting for simplicity E(q,p) — E and inverting the argument of the logarithm 



e- 0E S B/F (-q,-E) 



(2irh) 3 n(3q 1 - e~P E 
1 2irm 2 1 



log 



1 =F z exp 


\ ? ( 


-2mE+q 2 ) 2 ' 


8m 


q 2 


1 =F z exp 


( 


-2mE-q 2 ) 2 ' 


8m 


g 2 \ 



' {2irhf n(3q 1 



oflE 



log 



1 =F z exp 


/3 
8m 


(2mE+9 2 ) 2 " 
9 2 J 


1 =F z exp 



8m 


I 2 \ 



= S B/F (q,E), 



which proves (|12|). Similarly 



i-^Suni-^-E) 



1 2?rm 2 
(2tt^) 3 n/3g 

1 2?rm 2 
(2tt^) 3 n/3g 



z exp [— /3.E 1 ] exp 



/3 {-2mE + q 2 f 



■zexp 



8m 
[3 {2mE + q 2 ) 2 



8m 



S MB (q,E). 



The master equation (|l|) for the Rayleigh gas in the case of a free gas of Bose or Fermi particles takes therefore the 
form 



f t =M B/ .® 



= --[H ,q}+ / dn(q) 

ft JR 3 



(29) 



U(q)yJ S B/F {q,p)g^ S B/F (q,p)U ] (q) - -{S B/F (q,p),g} 



with S B/F (q,p) given explicitly by (p3[), and a similar expression A^ MB can be considered for a free gas of Maxwell- 
Boltzmann particles. Both M B/F and M MB are translation- and rotation-covariant and admit the same stationary 
solution with the canonical structure po(p). 



IV. BROWNIAN LIMIT 



In the framework of an ideal gas considered in Section [II, i.e., referring to M b /f and Ai M s, we now want to consider 



the physically distinguished case in which the test particle is much heavier than the particles making up the gas, so 
that a = m/M is much smaller than one, the so-called Brownian limit. In order to do this we have to evaluate the 
dynamic structure factor for a free gas in the case q<1. The natural starting points are expressions (27) and ( |2q ) 
in which the factor a has been put into evidence through the function (|2q). In particular we have the relations 



P ,2 



2m 



p-(q,p) 



—q 2 

8m 

P 2 

8m q 



JL 

2m 



(v(q,p)-q) 2 



2 2M 

2 
P 1 



_P_ 2 

8m 9 2 2M 



k' 2 

E(q,p) 

[q 2 



&-!*** 



= (3E(q,p). 



1 P 1 1 r 2 



2p ■ qf 



P™ p 2 

2q 2 



E 2 (q,p) 



i 1 1 r 2 ~ i2 

2p - q]+ 2f a ^ [q +2P ' q] 



P_™-p2 
2q 2 



E 2 {q,p) 



(30) 



The Brownian limit can now be taken neglecting the terms of order a in (30) or equivalently considering small energy 
transfer, corresponding to a broader time scale, and keeping in (|3(j) only the terms linear in E, disregarding higher 
powers of the energy transfer. The resulting dynamic structure factors, denoted by an index oo, are given by 



Sb/Aqip) 



1 2nm 2 



1 



(2irhf n(3q 1 _^[^ r+ a j f.] 



log 



1 =F ze~ 



-2 _0t_£_ 



g 2 L 2M ~ M 



1 =F ze~ 



JL„2 i0 r i IP 

8m^ C 2 \-2M ' M 



(31) 



and 



S^ B (q,p) = 



1 



2nm 2 



_JL„ 2 _0_I_qf_ , £P] 
8m ^ g 2 L2Af "i" M J 



(2tt?i) 3 n/3q 
respectively. Considering the corresponding expressions in terms of E(q, p 



(32) 



S~(q,E(q,p))=T 



1 2itm 2 



1 



(2irh) 3 n(3q 1 - e^^.P) 



log 



1 =F -ze" 



e 2 



lzpze-^e+^^'P) 



(33) 



and 



S^(g,£?(q,p)) 



1 27TTO 2 

(2tt?i) 3 n/3g 



■ze 8 



-&«%-#*(«*) 



(34) 



one immediately sees that rotational invariance is preserved in this approximation. One can check that also the 
detailed balance condition is not spoiled, in fact from ( |33| ) one has 



and from (| 



-pEqac , 

"b/pv 



-E) 



1 



T 



{2Trhf nfiq l-e-P E 
1 27rm 2 1 



log 



1 =F ze~ 



„ 2 , j8 p 



-, 0_„2 _0. 

1 =F Ze 8m! e 2 



(2tt?i) 3 n/3g 1 



/3B 



log 



1 =F ze sm« e 2 J 



1 =F -ze" 



„2 



? 2 e +fB 



= S^{q,E) 



e-^S~ B (-q,-E) 



1 27rm 2 



-ze s 



~£;<l 2 P -pE(q,p)&E(q,p 



(2ttH) 3 nf3q 
As a result in place of ( p9| ) we now consider the mapping Afcp 



e 2 ' 



^ B (<?,£)■ 



Tt =Ml 



M 



(35) 



-j:iHo,Q] 



dn(q) 



R ;i 



u{q)\l s^M,P)q\J s^MiP)^ {q) - 2 {5T/p(g»p)> q) 



ans similarly -Mj^ B for Maxwell-Boltzmann statistics. A^~ /F and -M^ B are still translation- and rotation-covariant 
and admit the same stationary solution with the canonical structure po(p). 

In the master equation (pi), or according to the physical system under consideration ( p9| ) or (pa), the quantum 
scattering rate or transition probability appears through the dynamic structure factor and the square modulus of 
the Fourier transform of the T matrix determining the integration measure (0), these quantities being connected 
to the differential scattering cross-section by (0). In order to pass from the master equation to the related Fokker- 
Planck equation through a Kramers-Moyal expansion, as stressed by van Kampen p7| we need to put into evidence 
a small parameter governing the size of the fluctuations in the macroscopic system. In our case this parameter is 
naturally given by the momentum transfer q, which through the dynamic structure factor is directly linked to the 
equilibrium fluctuations of the macroscopic system. Small q means long- wavelength fluctuations, corresponding to the 
macroscopic, long range properties of the environment. It is physically meaningful to consider both approximations 
\q\ <§C 1 and a«l, or equivalently small energy transfer, together, so that starting from the Maxwell-Boltzmann case 



i ~ Air tyi 

•MmbH = - jWo, ] + z—rr- I d 6 q -^ ^< 



3 \i(q)\ 2 j^ (1+2a)Q 2 



\ {.-*-,} 



we expand the dissipative part of the mapping in the small parameter q. We will expand the operators depending on 
q up to second order, so as to have contributions at most bilinear in the operators x and p, position and momentum of 
the Brownian particle. We thus obtain a structure analogous to the classical Fokker-Planck equation, with a friction 
term linearly proportional to velocity: this class of models is known as quantum Brownian motion |1S}-|2C|]. Recalling 
that a<l and exploiting the symmetry properties of the integration measure the result for the dissipative part is 17]] 



2tt 2 to 2 
0h 



^ML-^ 2 



-e 8" 



R3 



H^{^2 [«i. [«*)•]] 



I j^\p i ,\P i ,-]} + ~^l*i,{Pi,-}] 



where z = 1,2,3 denotes Cartesian coordinates. Due to isotropy of the environment we have q\ = -|(7 2 , so that we can 
define the coefficients 



D n 



D, 



_ 2 7T 2 TO 2 

~ 3 j3h 
\2M 



R3 



d 3 q\i(q)\ 2 qe- 



S- a 2 
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(36) 
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and introduce the following mapping describing quantum dissipation 



A,d[-1 = - : 



D 



Ei 



D. 
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(37) 



i=l 



thus coming to the Fokker-Planck equation 



do 



-[H ,g] +C MB [g] 



k 



[Ho,e] + zC QD [g\. 



(38) 



In view of the result ( pq ) for the Fokker-Planck equation describing the motion of the Brownian particle in a gas 
obeying Maxwell-Boltzmann statistics, we now look for the corrections to ( p8[ ) brought about by quantum statistics 
at finite temperature. As usual we will deal with both Bose and Fermi statistics, exploiting expression (Bq) obtained 
in Appendix [b| by deriving an exact expansion for S™ /F (q, E): 



S-( q> E) = S^(q,E) 



E^-^^Ew^ 



Lfc=0 



(39) 



fc=l 



OC 1 OO 



fc=i 



fc=i 



where a suitable expansion in the small parameter q has already been performed. Let us first introduce the Bose- 
Einstein and the Fermi-Dirac functions pi 22 1, given by 



and 



9u{z) 



fu{z) = 



1 



T{v) 



dx ■ 



z 1 e x — 1 



-K— 1 



dx ■ 



T{u) J z- 1 e x + 1 



0<z<l,u>0 



< z < oo, ^ > 



(40) 



(41) 



respectively, related for integer z; by f n (z) = —9n{—z). These functions, typically appearing in the quantum statistical 
mechanics of Bose and Fermi systems, satisfy the recurrence relations 



d d 

9u-i{z) = z—[g v (z)\ f v -x{z) = z— \f v {z) , 
az oz 



(42) 



so that they can be defined also for v < 0. Starting from (E2|) and exploiting the following representations for \z\ < 1 



00 u 
z K 



fj, 



w = Et /^) = E(-) fe - 1 fc. 



A,-i 



fc=i 



one can write (|39j) in the Bose case as 
S~(q,E) = S~ B (q,E) 

and in the Fermi case as 

S~(q,E) = S~ B (q,E) 



""'' : ' 3 q 2 9-i(z) + i-MEf(9-2{z) + 2g_ x {z)) + 0( 9 4 ) 



8m 



1 
24" 



' / " ( ' ;) + ^-q 2 .f-i(z) - i-([3E) 2 {f-2{z) + V-i{z)) + 0(<z 4 ) 



1 

24' 



Z 8m 

where the functions appearing in ( [43[ ) can be written for < z < 1 in closed form as: 

9o(z) 

5 - l(z)= (I^ 5 - 2(z)= (l-z)3' 
while the functions appearing in (E3) can be written for < z < 00 in closed form as: 

/o(*) = L 



z 
\-z 

z 



(43) 



(44) 



(45) 
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1 + z 
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(46) 
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To evaluate the corrections due to quantum statistics we note that when S^ B (q, E{q, p)) is substituted by an expression 
of the form 

Smb(<7> E(q,p))A[l + 2Bq 2 + 2C(p ■ q)% 

so that keeping terms at most quadratic in q in the correction 

V5» B (g,S(g,p)) - v/5» B (g,£;(g,p))VI[l + Bg 2 + C(p • g) 2 ], 
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the mapping -Mj^ B always in the same approximation becomes 



tfHo,.] 



f A I dn(q) 

'R3 



(47) 
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U(qWS% B (q,P){(p • q)\ -}y/ ' S£ B {q,P)U\q) - {5» B (g,p)(p • g) 2 , •} 



Looking at (47) one immediately sees that, keeping terms at most quadratic in q, the last two terms are to be neglected, 
since the dynamic structure factor S^ B (q,p) and the unitary operators U(q) can now only bring in a constant factor. 
The only change in the structure of the mapping is therefore given by the numerical factor A multiplying the dissipative 
part. This factor is actually given by g (z)/z in the Bose case and by f (z)/z in the Fermi case. The Fokker-Planck 
equation (Bq) in the case of Bose statistics of the gas therefore becomes 



dp 
~dl 

while for Fermi particles one has 



■[H , g] + C b [q] 



d j t =^-[H Qie] + c M 

and the following simple relations hold 

£b = 5o( z )-£qd 



-[H ,q] +go(z)£ Qn [g] 



-r[-Ho,e] + fo(z)C QD [g], 



(48) 
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According to (pfl) and setting after (Bq) 



/3 1 7r 2 m 2 



1 + 2 



-73„ 1 + ^12^-^q 
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one has the following very simple relation between the friction coefficients in ( pq) and 

Tmb 



or (49) 



7b/f 



lT«' 



(51) 



(52) 



The relationship between the Fokker-Planck equations for Maxwell-Boltzmann or Bose and Fermi statistics, as given 
respectively by (pq), ( ^8[ ) and (p9|), is actually remarkably simple: they have the very same operator structure, apart 
from an overall coefficient depending on the fugacity of the gas, which determines the relative weight of the dissipative 
contribution to the dynamics. As it is to be expected only the statistics of the reservoir is of relevance, since the test 
particle is a single particle. The Fokker-Planck equations obtained for the description of quantum dissipation may be 
compactly written: 



do 
~dl 



-[H ,q] + ((z)C q 4q] 



with £(z) defined as follows 



C(z) 



z Maxwell — Boltzmann 

z/(l — z) Bose 

z/(l + z) Fermi 



(53) 



(54) 
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We now briefly come back to the question dealt with at the end of Section [Fj] about the physical relevance of the 
structure of translation-covariant master equation obtained in [|16| . As already stressed the master equation (|l|) , while 
having t he g eneral translation-covariant structure considered in pl[ , does not comply with the further restrictions 
given in fLq] , while this would be the case if instead of the dynamic structure factor S(q,p) one would consider the 
symmetrized correlation function S(q,p), which is an even function of q and E(q,p). This could be considered a 
natural phenomenological Ansatz in view of the result (l^] . In Section |l] we showed however that this substitution 
would spoil the existence of the expected stationary solution. More than this, if we now consider the Brownian limit, 
the symmetrized version of S^ B (q,p), which can be immediately obtained from (|34), reads 



S~ B (q,E(q,p)) = 



1 2nm 2 
(2Trh) 3 n[3q 



so that the dependence on p is completely lost and the whole operator structure in (|37]) and ( |38| ) is washed out, apart 
from the double commutator with the position operators Xi. In particular the friction term is missing, so that, even 
though a Lindblad structure is retained, only a completely different physics can be described. In the same way, for 
the Bose or Fermi dynamic structure factor in the Brownian limit one has from (B6) 



S-(q,E)=S- B (q,E) 



8tti ^ k ' 

k=l 



So**) 

fc=0 

-i(3E) 2 E(±) fcfc ^ + ^ E ) 2 E(±) fcfcV + °(« 4 ) 
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k=l 



24 



k=i 



and once again, recalling ( |47| ) written in terms of S^ B rather than S^ B , one sees that under the same approximations 
as before the operator structure in the dissipative part of both (ES) and (Bfl) is washed out apart from the contribution 
due to the double commutator in the position operators of the particle Xi . 

We now consider some structural features of the mapping C QD given by (pTj) in terms of which the Fokker-Planck 
equation ( p3[ ) encompassing all three statistics is given. G-covariance of £ QD under translations and rotations imme- 
diately follows from its very structure and the transformation laws for the operators x and p: 

U ] {a)xU{a) = x + a, U f (a)pU(a) = p, U\R)xU{R) = Rx, U\R)pU(R) = Rp. 

One can also see that an operator with the expected canonical structure is a stationary solution of (|52|) in that 

£qdMp)] = o, 

due to the relationship 
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2M 



(55) 



obeyed by the coefficients defined in ( |3q ) and entering in (37). A few more remarks are in order. The typical structure 
of translation-covariant mappings describing quantum dissipation in analogy with the classical Fokker-Planck equation 
that one finds in the physical literature is given by p3| 
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(56) 



where the ratio between the first two coefficients, given by (3/2M as in (|5^) is fixed by the requirement that po(p) be 
a stationary solution, i.e., £p P [po(p)] == 0, and the only freedom left, apart from the overall multiplying coefficient 7 
is given by the adimensional factor \. If one further asks that ([36]) can be cast in Lindblad form, so that £p P is the 
generator of a completely positive dynamical semigroup Q , one has the further simple requirement p4],|l8[ 



X 



1 
> -. 



(57) 



In fact under this condition, observing that for the operators 

Bi± = Xi ± inpi 
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we have the identity 



- t 1 " t - 1 

B l± B t± - -{B l ±B l ± 1 -} = —z{[xi,[xi,'}] + k 2 [p^ipi, ■}} ± 2iK[x l7 {p l7 ■}} T in[{x l ,p i } ,•]}, 



we may write £p P in an explicit Lindblad form in terms of the two generators 

f . ,M fx. t 

L i+ = Xi + i—\l —p.-, Li- = x 
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(58) 



The Fokker-Planck structure £ QD falls within this class, with the coefficient 7 given by ( pq) in terms of a suitable 
integral of the Fourier transform of the T matrix describing collisions at microphysical level. Moreover it corresponds 
to the value \ = 1/8 in (|57|), so that 



r _ -1/8 

*-QD — <-FP ■ 



(59) 



This in turn implies that £ QD can be written in a manifest Lindblad form in terms of a single generator for each 
Cartesian direction. We make the choice ftpi! 



V2 /„ i A 2 



A M V X< ' fi 4 



M, 



where Am = \ h (3/M, the thermal wavelength associated to the Brownian particle, is put into evidence, so that one 
has the commutation relations 



In such a way we have the alternative expression 

12 3 



^o[-] = -^^£J[{^pJ,-] + ^a m ]T 



i=l 



in which the single generator feature is put into evidence. 



Oi-al- -{a\a,i,-} 



(60) 



V. CONCLUSIONS AND REMARKS 



In this paper we have considered the behavior with respect to covariance under translations and rotations, and 
existence of a stationary solution of a recently proposed master equation (|l|) for the description of the interaction of 
a test particle with a fluid, a physical model corresponding to the so-called Rayleigh gas. The key result in (Q) is the 
appearance of a two-point correlation function known as dynamic structure factor and given by (^) , the general struc- 
ture conforming to results already obtained in the mathematical literature for the generator of a translation-covariant 
dynamical semigroup. This correlation function, depending on symmetry and statistical mechanics properties of the 
fluid, directly determines the behavior of the master equation with respect to covariance under translations and rota- 
tions, and existence of a stationary solution with the expected canonical form. Considering the specific case of a free 
gas, the dynamic structure factor has been explicitly evaluated for Bose, Fermi and Maxwell-Boltzmann statistics, and 
the dependence on the physical parameters determining the peculiar features of the model under consideration has 
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been put into evidence in an exact expansion of the dynamic structure function. These parameters are the fugacity 
of the gas z, the ratio between mass of the gas particles and of the test particle en, the transferred momentum q and 
the transferred energy E(q,p). Stability of the covariance properties of the master equation and of existence of a 
stationary solution is then considered in the limit in which these parameters are small, together with the different 
explicit expressions of the master equation. In particular in the Brownian limit a«l and considering small momen- 
tum transfer, corresponding through the physical interpretation of the dynamic structure factor to long-wavelength 
fluctuations, a Fokker-Planck equation with a Lindblad structure is obtained, given by (p3[), where the results cor- 
responding to Bose, Fermi and Maxwell-Boltzmann statistics are jointly considered. The correction due to quantum 
statistics in the Fokker-Planck equation is simply expressed through the Bose and Fermi functions given by ( fig ) and 
(|46|) respectively. 
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APPENDIX A: DERIVATION OF EQ. (E3) AND EQ. (25) 



In this Appendix we want to explicitly calculate the expression of the dynamic structure factor for a free gas as a 
function of q and p. Working at finite temperature we can carry out the calculation for both Bose and Fermi particles 
at the same time adopting the convention that the symbol ± means a + sign for Bose particles and — for Fermi 
particles. We start from the expression 

<? i nn \ 1 f rf3fc / \ m + / \ ^ (p + <?) 2 , ( fc -q) 2 p 2 k2 \ 



with 



1 k 2 

z L e pek =F 1 2m 



which can be found for example in p2[ and corresponds to (0) for a free gas apart from a singular term proportional, 
in the continuum limit considered here, to S 3 (q), relevant only for q = and not contributing to the master equation. 
In fact in the derivation of the master equation the contributions for q = exactly cancel out. This term according 
to (fn) corresponds to forward scattering. We now have to evaluate the integral in k. This is most easily done writing 
Sb/f(q,p) in the form 



and observing that 



1 f d k 
Sn/v{q,P) = - I ,„ rrg (n k ) B/F (l ± {n k ^ q ) B/F )5 (E(q,p) + e k - 
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so that one has 
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We are thus led to consider an integral of the form 

/ d 3 k(n k ) B/F 5(ri-2k- q) 

JR3 



(A2) 



with 77 a real parameter. Denoting by £ the cosine of the angle between k and q the integral in ( A2) becomes 
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and using the identity 
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we get, integrating by parts, 
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Inserting the result ( |A5| ) in (Al) one immediately has 
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In a similar way starting from the expression of a gas of Maxwell-Boltzmann particles 
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we write it in the form 
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Analogously to (A2) we have to consider 



d 3 k(n k ) MB 5(i]- 2k ■ q) , 



which according to (A3) becomes 
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Exploiting instead of (A4) the relation 
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i / \ m d -e« 
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we obtain 
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which has to be substituted in (A7) leading to 
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APPENDIX B: EXACT EXPANSION OF Sb/f AND DERIVATION OF EQ. (B9) 



We will now derive an expression for S B/F (q,p) equivalent to (|2j]), in which however a series expansion in powers 
of the fugacity z is put into evidence. The starting point is the Taylor expansion for the logarithm log(l + x) — 
J^kLi {—) k+1 x k /k, which leads to write ( p7| ) in the form 
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Considering now in (Bl) a geometrical progression of reason e2>™( 2<T (<z<ph 9 ', according to the formula 



fc-i 



1-x* 



(i -*)$>" 



(B2) 



n=0 



cq. (Bl) becomes 
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where it is to be noted that the sum over n is due to the presence of the statistical correction 1 ± (nk- q ) B / F hi ( |22| ) 
and disappears, being replaced by a factor one, if this correction is neglected. It is worthwhile to put into evidence a 
factor 



1 27rm 2 
(2tt?i) 3 n/3q '' 



-£:- 2 (i-p) 



corresponding to the expression of the dynamic structure factor for a free gas of Maxwell-Boltzmann particles, thus 
obtaining 



S B/F (q,p) = S MB (q,p) 
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In the Brownian limit a< 1 considered in Section IV, neglecting in ( p0[ ) the contributions of order a, this expression 
goes simply over to 



S B/F (q,p) = S~ B (q,p) 






fe=i 
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(B3) 



where S B/F (q,p) and S^ B (q,p) are given respectively by ( |3l"| ) and (|32"|). In terms of E(q,p) equation (B3) takes the 
remarkably compact form 
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where S B/F (q,E) and S^ B (q,E) are given by (p3|) and (p4). We now go one step further noting that the following 



identity holds: 
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which can be easily obtained exploiting (B2), so that (B4) becomes 
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Eq. (B5) is the most convenient expression in order to consider the limit of small momentum transfer. Exploiting the 
expansion 



sinh(±/3£) V ; 
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and recalling that E is given by (|6|) we may write (B5) as 



S~(q,E) = SZ B (q,E) 
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Recalling (|B5| ) and the explicit expression of 5~ B given by (£2|) one also has 
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we obtain the alternative expression 
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equivalent to (|3^) as can also be directly checked starting from the identity 
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Note that ( pl|) and (B7) in the Fermi case can also be written in the form 
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and 
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^ / ^n 1 27TTO 2 < - ~ 

S^(q,E) = - . ... — ? r arth 

(2^)3 n/3g sinh (|^ 

respectively, that can be useful if one is interested in an expansion for large values of z. According to (BE) also for 
S B/F (q,E) given by (|23| ) one has the alternative expression 



^-fl-e-^^coshUE)) I 



1 27rm 2 e 2 B 
arth 



i mm e = 
{2-Khf n/3q sinh ^ E 



±ze-s^9 2 e " 2 ^ E2 smh(§E 



-^„-^ E2 



1 =F ze »m« e 



cosh f 



l£ 



(B9) 



balance condition for ( p9| ) according to (fi"^) can immediately be checked observing that 
eld functions, while coshx is an even function. 



The validity of the detailed uaiaiitu ^uiiuiliuh iui \v^\j c 
both sinh a: and arthx are odd functions, while cosh a: is an 
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